6 years of harvest, 2020 had variable rate seeding

2013

plot.gps.grids(corn.2013.dat,plots=0,residuals=TRUE) 
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.

2014

plot.gps.grids(soybean.2014.dat,plots=0,residuals=TRUE) 

2015

plot.gps.grids(corn.2015.dat,plots=0,residuals=TRUE) 

2016

plot.gps.grids(soybean.2016.dat,plots=0,residuals=TRUE) 

2017

plot.gps.grids(corn.2017.dat,plots=0,residuals=TRUE) 

2018

plot.gps.grids(soybean.2018.dat,plots=0,residuals=TRUE) 

Variograms

Yield Monitor

Full Field

corn.2013.var <- variogram(Yield~1, 
                          locations=~Easting+Northing,
                          alpha=c(0,45,90),
                          data=corn.2013.dat)
plot(corn.2013.var)

soybean.2014.var <- variogram(Yield~1, 
                          locations=~Easting+Northing,
                          alpha=c(0,45,90),
                          data=soybean.2014.dat)
plot(soybean.2014.var)

corn.2015.var <- variogram(Yield~1, 
                          locations=~Easting+Northing,
                          alpha=c(0,45,90),
                          data=corn.2015.dat)
plot(corn.2015.var)

soybean.2016.var <- variogram(Yield~1, 
                          locations=~Easting+Northing,
                          alpha=c(0,45,90),
                          data=soybean.2016.dat)
plot(soybean.2016.var)

corn.2017.var <- variogram(Yield~1, 
                          locations=~Easting+Northing,
                          alpha=c(0,45,90),
                          data=corn.2017.dat)
plot(corn.2017.var)

soybean.2018.var <- variogram(Yield~1, 
                          locations=~Easting+Northing,
                          alpha=c(0,45,90),
                          data=soybean.2018.dat)
plot(soybean.2018.var)

References